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Summary. We review a recent global analysis of inclusive single-charged-hadron 
production in high-energy colliding-beam experiments, which is performed at next- 
to-leading order (NLO) in the parton model of quantum chromodynamics endowed 
with nonperturbative fragmentation functions (FFs). Comparisons of pp data from 
CERN SppS and the Fermilab Tevatron and yp data from DESY HERA with the 
corresponding NLO predictions allow for quantitative tests of the scaling violations 
in the FFs and their universality. We emphasize the potential of new measurements 
at the Tevatron to place tight constraints in the large-a; region and on the gluon 
FF, complementary to those from e + e _ data, which are indispensible in order to 
reliably predict the 7 r° background for the H —> 77 signal of the intermediate-mass 
Higgs boson at the Tevatron and the CERN LHC. Adopting a similar theoreti¬ 
cal framework for 6 -hadron production, we show that the notorious excess of the 
Tevatron data over existing theoretical calculations can be ascribed, at sufficiently 
large values of pr, to nonperturbative fragmentation effects inadequately included 
previously. 


1 Introduction 

In the framework of the QCD-improved parton model, the inclusive produc¬ 
tion of single hadrons is described by means of fragmentation functions (FFs) 
D%(x, p 2 ). The value of D%(x, p 2 ) corresponds to the probability for the par- 
ton a produced at short distance 1/p to form a jet that includes the hadron 
h carrying the fraction x of the longitudinal momentum of a. Unfortunately, 
it is not yet possible to calculate the FFs from first principles, in particular 
for hadrons with masses smaller than or comparable to the asymptotic scale 
parameter A. However, given their x dependence at some energy scale p, the 
evolution with p may be computed perturbatively in QCD using the timelike 
Altarelli-Parisi (AP) equations |jj. This allows us to test QCD quantitatively 
within one experiment observing single hadrons at different values of centre- 
of-mass (CM) energy yfs (in the case of e + e _ annihilation) or transverse 
momentum pp (in the case of scattering). Moreover, the factorization theo¬ 
rem guarantees that the D^(x,p 2 ) functions are independent of the process 
in which they have been determined and represent a universal property of h. 

* To appear in the Proceedings of the lfth Topical Conference on Hadron Collider 
Physics (HCP 2002), 29 September - 4 October 2002, Karlsruhe, Germany. 
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This enables us to make quantitative predictions for other types of experi¬ 
ments as well. 

During the last seven years, the experiments at CERN LEP1 and SLAC 
SLC have provided us with a wealth of high-precision information on how par- 
tons fragment into low-mass charged hadrons (h^). The data partly comes 
as light-, c-, and 6-quark-enriched samples without or with identified 
final-state hadrons (tt^ , K ± , and p/p) [Q or as gluon-tagged three-jet sam¬ 
ples without hadron identification |5|,|(J. Motivated by this new situation, the 
author, together with Kramer and Potter, recently updated, refined, and 
extended a previous analysis (t| by generating new leading-order (LO) and 
next-to-leading-order (NLO) sets of 7r ± , I\ ± 1 and p/p FFs By also in¬ 
cluding in our fits 7r , A , and p/p data (without flavour separation) from 
PEP [jjj, with CM energy yfs = 29 GeV, we obtained a handle on the scal¬ 
ing violations in the FFs, which enabled us to determine the strong-coupling 
constant o$\Mz )■ 

The formation of D and B mesons from c and 6 quarks, respectively, 
is a genuinely nonperturbative process, so that it is indispensable to em¬ 
ploy nonperturbative FFs. Furthermore, if the characteristic energy scale p 
is large against the heavy-quark mass m c ( mb ), then the QCD-improved par- 
ton model with nf = 4 (n/ = 5) active quark flavours is the appropriate 
theoretical framework, which allows for coherent analyses of data from e + e _ , 
lepton-hadron, and hadron-hadron collisions | fTo||TT) . 

This contribution is organized as follows. In Sect. we present some 
details of our global fits (8) and assess the quality of the resulting FFs. We 
also discuss the determination of ai 5 \Mz) from the scaling violations in 
the FFs |0]. In Sect. I we present comparisons of our NLO predictions for 
inclusive charged-hadron production |k| with pp data from SppS 0 and the 
Tevatron [[[ 5 ] and with yp data from HERA |I(|. In Sect. [|, we demonstrate 
that the pt distribution of inclusive R-meson hadroproduction measured by 
CDF 0,^8) is well described at NLO in the QCD parton model with nf = 5 
endowed with nonperturbative FFs, once the latter are fitted to LEP1 data 
[p~9| . Our conclusions are summarized in Sect. ||. 


2 Determination of the FFs 

The NLO formalism for extracting FFs from measurements of the cross sec¬ 
tion da/dx of e + e _ —> h + X is comprehensively described in [j7j. We work 
in the modified minimal-subtraction (MS) renormalization and factoriza¬ 
tion scheme and choose the renormalization and factorization scales to be 
Pr = Pf = £\/s, except for gluon-tagged three-jet events, where we put 
PR = Pf = £ x 2£j e t, with £jet being the gluon-jet energy in the CM frame. 
Here, the dimensionless parameter £ is introduced to determine the theoret¬ 
ical uncertainty in ai 5 \Mz) from scale variations. As usual, we allow for 
variations of £ between 1/2 and 2 about the default value 1. 
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Our strategy was to only include in our fits LEP1 and SLC data with 
both flavour separation and hadron identification [Q, gluon-tagged three- 
jet samples with a fixed gluon-jet energy |f|, and the 7r ± , A' ± , and p/p 
data sets from the pre-LEPl/SLC era with the highest statistics and the 
finest binning in x j9|. The y 2 value per fitted data point turned out to be 
Xdf = 0.97 (0.98) at LO (NLO). The goodness of our fit may also be judged 
from Figs. |l](a) and (b), where our LO and NLO fit results are compared 
with the ALEPH, DELPHI, OPAL, and SLD data • Other data served 
us for cross checks [||j|,|],|o],^l] . In particular, we probed the scaling viola¬ 
tions in the FFs through comparisons with 7r , A , and p/p data from DESY 
DORIS and PETRA, with CM energies between 5.4 and 34 GeV pof . Fur¬ 
thermore, we tested the gluon FF, which enters the unpolarized cross section 
only at NLO, by comparing our predictions for the longitudinal cross section, 
where it already enters at LO, with available data Finally, we directly 

compared our gluon FF with the one recently measured by DELPHI in three- 
jet production with gluon identification as a function of x at various energy 
scales p [(|. All these comparisons led to rather encouraging results. We also 
verified that our FFs satisfy reasonably well the momentum sum rules, which 
we did not impose as constraints on our fits. 

Since we included in our fits high-quality data from two very different en¬ 
ergies, namely 29 and 91.2 GeV, we were sensitive to the running of a s (p) and, 
therefore, able to extract values of A—. We obtained A^—L = 88 f + , 3 o MeV 

at LO and A^—L = 213 + 1\ + ^ 2 MeV at NLO, where the first errors are exper- 
imental and the second ones are theoretical. From the LO and NLO formulas 
for we thus obtained af\Mz) = 0.1181 ^ooo 69 -o 0049 (LO) and 

ai 5 \Mz) = 0.1170^o Qggg I 0 0025 (NLO), respectively, which is already in¬ 
cluded in the most recent world average by the Particle Data Group p2{ . 
As expected, the theoretical error on a., 5) (Mz) is significantly reduced as we 
pass from LO to NLO. We observe that our LO and NLO values of ai 5 \Mz) 
are quite consistent with each other, which indicates that our analysis is per- 
turbatively stable. The fact that the respective values of ALA. significantly 

differ is a well-known feature of the MS definition of a[ nf \p) [f23| . 


3 Global Analysis of Collider Data 

Recently, we extended our previous tests of scaling violations (|| to higher 
energy scales by confronting new data of e + e _ —^h ± +X from LEP2 J24|, 
with y/s ranging from 133 GeV up to 189 GeV, with NLO predictions based 
on our FFs Jl3| . Furthermore, we quantitatively checked the universality of 
our FFs by making comparisons with essentially all available high-statistics 
data on inclusive charged-hadron production in colliding-beam experiments 
[ p~3| . This includes pp data from the UAl and UA2 Collaborations |L4| at 
SppS and from the CDF Collaboration [li| at the Tevatron, yp data from 
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Fig. 1. (a) Comparison of data on inclusive charged-hadron production at y/s = 
91.2 GeV from ALEPH ( triangles ), DELPHI (circles), and SLD ( squares ) |Ij with 
our LO (dashed lines) and NLO (solid lines) fit results 0. The upmost, second, 
third, and lowest curves refer to charged hadrons, 7r , A , and p/p, respectively, (b) 
Comparison of three-jet data on the gluon FF from ALEPH with Ej e t = 26.2 GeV 
(upper curves) and from OPAL with Ej e t = 40.1 GeV (lower curves) 0 with our 
LO (dashed lines) and NLO (solid lines) fit results 0 


the HI and ZEUS Collaborations JuJ at HERA, and 77 data from the OPAL 
Collaboration ^5| at LEP2. In hadroproduction and photoproduction, we set 
[i r = /if = £pt- As for the parton density functions (PDFs) of the proton, 
we employed set CTEQ5M provided by the CTEQ Collaboration |26|, with 
A— = 226 MeV. As for the photon PDFs, we used the set by Aurenche, 
Fontannaz, and Guillet (AFG) [E7|. In all cases, we found reasonable agree¬ 
ment between the experimental data and our NLO predictions as for both 
normalization and shape, as may be seen from Fig. [| We conclude that 
our global analysis of inclusive charged-hadron production provides evidence 
that both the predicted scaling violations and the universality of the FFs are 
realized in nature. 

4 Inclusive 13-Meson Production 

The QCD-improved parton model implemented in the MS renormalization 
and factorization scheme and endowed with nonperturbative FFs, which 
proved itself so convincingly for light- and D* ± -meson [|lo) inclu¬ 

sive production, also provides an ideal theoretical framework for a coherent 
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Fig. 2. Comparisons of (a, b) SppS |l4| and (c) Tevatron jTbJ data of pp —» ft* + X 
and (d) HERA data of 'yp —> h ± + X |1g|] with our NLO predictions 


global analysis of .B-meson data ||Tl[, provided that /i> mb, where ji is the 
energy scale characteristic for the respective production process. Then, at 
LO (NLO), the dominant logarithmic terms, of the form a"’" +1 In" (/r 2 /m^) 
with n = 1,2,..., are properly resummed to all orders by the AP evolution, 
while power terms of the form (m 2 /p. 2 ) are negligibly small and can be 
safely neglected. The criterion p, 3> mb is certainly satisfied for e + e _ an¬ 
nihilation on the Z -boson resonance, and for hadroproduction of B mesons 
with pp mb- Furthermore, the universality of the FFs is guaranteed by 
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the factorization theorem |2§|] , which entitles us to transfer information on 
how b quarks hadronize to B mesons in a well-defined quantitative way from 
e + e _ annihilation, where the measurements are usually most precise ]l9|j29| |, 
to other kinds of experiments, such as hadroproduction Ell¬ 
in |ijJ , the distribution in the scaled B -meson energy x = 2 Eb / y/s mea¬ 
sured by OPAL |l^] at LEP1, which is compatible with the subsequent mea¬ 
surements by ALEPH at LEP1 and SLD at SLC ]29| , was fitted at LO and 
NLO using three different ansatze for the b —> B FF at the starting scale 
p 0 = 2mb = 10 GeV. The best fit was obtained for the ansatz by Peterson 
et al. (P) Qi with Xdf = 0.67 (0.27) at LO (NLO) [see Fig. |](a)]. The e 
parameter was found to be e = 0.0126 (0.0198). We emphasize that the value 
of e carries no meaning by itself, but it depends on the underlying theory 
for the description of the fragmentation process b —> B 1 in particular, on 
the choice of the starting scale po, on whether the analysis is performed in 
LO or NLO, and on how the final-state collinear singularities are factorized 
in NLO. In Fig. |](b), the B + - m eson pr distribution measured by CDF in 
Run 1A jlTj] and Run 1 |l^| is compared with our LO and NLO predictions 
evaluated using the CTEQ6 proton PDFs El and the BKK-P B-meson FFs 
[ IN] and choosing pr.f = £r,f x 2 rriT, where tot = \Jl?T + m b- The the¬ 
oretical uncertainty at NLO is estimated by independently varying £r and 
from 0.5 to 2 about the default value 1. It only amounts to at 

Pt = 15 GeV, but steadily increases towards smaller values of pt, reaching 
(*( 4 g% at pt = 10 GeV. Variations in the proton PDFs and the ansatz for the 
fc^BFF 0 only reach a few percent. We observe that the Run 1A data 
comfortably lies within the theoretical error band, while the full Run 1 data 
tends to be somewhat on the high side, especially in the upmost pt bin. 

In the case of 77 —> D* ± +X at LEP2, the inclusion of finite-m c effects was 
found to reduce the cross section by approximately 20% (10%) at pr = 2 m c 
(3 m c ) E|, i.e., their magnitude is roughly as naively expected. By 

analogy, one expects the finite-7?i6 terms neglected in JdJ to have a moderate 
size, of order 20% (10%) at pr = 10 GeV (15 GeV). This is considerably 
smaller than the scale uncertainty and appears insignificant compared to the 
excess of the CDF data [FtUr] over the traditional NLO analysis |k| in a 
scheme with n/ = 4 active quark flavours, where the b quark is treated in 
the on-mass-shell renormalization scheme. In this massive scheme, there are 
no collinear singularities associated with the outgoing 6-quark lines that need 
to be subtracted and absorbed into the FFs. In fact, in this scheme there is 
no conceptual necessity for FFs at all, but they are nevertheless introduced 
in an ad-hoc fashion in an attempt to match the B-meson data. However, in 
the absence of a subtraction procedure, there is also no factorization theorem 
in operation to guarantee the universality of the FFs E|. By the same to¬ 
ken, such FFs are not subject to AP evolution. Thus, there is no theoretical 
justification to expect, e.g., that a single value of the Peterson e parameter 
should be appropriate for different types of experiment or at different energy 
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(b) 



Pt [GeV] 


Fig. 3. (a) Our LO and NLO fits @ to OPAL data of e+e" —> B + X (b) 
Comparison of CDF data of pp —> B + X E3JH] with our LO and NLO predictions 

0 


scales in the same type of experiment. In other words, the feasibility of global 
data analyses is questionable in this scheme. Moreover, this scheme breaks 
down for px m& because of would-be collinear singularities of the form 
a s In , which are not resummed. 

The attempt to split the B-meson FFs into a so-called perturbative FF 
(PFF) and a nonperturbative remainder is interesting in its own right. How¬ 
ever, detailed analysis for D* ± -meson FFs Q revealed that such a procedure 
leads to deficient results in practical applications. On the one hand, at NLO, 
the cross section da/dx of e + e _ annihilation becomes negative in the upper 
x range, at x >0.9 [see Fig. [|(a)], where the data is very precise, so that a 
low-quality ht is obtained unless this x range is excluded by hand . On 

the other hand, the LO and NLO predictions for other types of processes, 
such as photoproduction in ep scattering [see Fig. [|(b)], significantly differ 
[0| , which implies that the perturbative stability is insufficient. 

The idea jnj of performing a coherent analysis of LEP1 and Tevatron 
data of inclusive B-meson production was recently revived using an uncon¬ 
ventional scheme named FONLL, in which the traditional result in the mas¬ 
sive scheme [0J and a suitably subtracted result in a massless scheme with 
PFFs are linearly combined |3(j]. The degree of arbitrariness in this procedure 
may be assessed by noticing that the massless-scheme term is weighted with 
an ad-hoc coefficient function of the form p\/ (pj, + 25 mf) so as to effectu¬ 
ate its suppression in the low-p-r range and that this term is evaluated at 
p' T = \/p\ + while the massive-scheme term is evaluated at px- Since the 
FONLL scheme interpolates between the massive scheme and the massless 
scheme with PFFs, it inherits the weaknesses of both schemes detailed above. 
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Fig. 4. (a) NLO fit §| to OPAL data of e+e~ -*■ D*± + A' using PFFs. (b) 
Comparison of ZEUS data of ep —> D* ± + X in photoproduction with our LO and 
NLO predictions implemented with PFFs Q 


In particular, the negativity of the NLO cross section of e + e~ —> B + X in 
the upper x range forces one to exclude the data points located there from 
the fit. In |l6|, this is achieved by resorting to what is called there the mo¬ 
ments method, i.e., the large-cc region is manually faded out by selecting one 
particular low moment of the b —> B FF, namely the one corresponding to 
the average x value, thereby leaving the residual information encoded in the 
data unused. 

5 Conclusions 

We reviewed recent LO and NLO analyses of 7r ± , K ± , and p/p FFs ||, which 
also yielded new values for a^\M z ) ||l2). Although these FFs are genuinely 
nonperturbative objects, they possess two important properties that follow 
from perturbative considerations within the QCD-improved parton model 
and are amenable to experimental tests, namely scaling violations and uni¬ 
versality. The scaling violations were tested |^,^3| by making comparisons 
with data of e + e _ annihilation at CM energies below |^o) and above [^4j 
those pertaining to the data that entered the fits. The universality prop¬ 
erty was checked ||l3) by performing a global study of high-energy data on 
hadroproduction in pp collisions [H| . [i~T)| and on photoproduction in e ± p [ pu| 
and e + e _ (25) collisions. Our NLO FFs Q agree with other up-to-date sets 
| [l7| within the present experimental errors |p| . 

High-energy data on hadroproduction particularly probes the FFs in the 
large-cc region, is especially sensitive to the gluon FFs, and, therefore, carries 
valuable information complementary to that provided by e + e~ data. This 
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information, apart from being interesting in its own right, is indispensible 
in order to reliably predict the 7 r° background for the H —> 77 signal of 
the intermediate-mass Higgs boson at the Tevatron and the LHC. It would 
thus be highly desirable if the experiments at the Tevatron made an effort to 
update their analyses, which date back to 1988 Jl5| ]. 

Adopting a similar theoretical framework for H-meson production, we 
demonstrated [jll) that the notorious excess of the Tevatron data 0 ® over 
existing theoretical calculations can be ascribed, at sufficiently large values 
of pt , to nonperturbative fragmentation effects inappropriately taken into 
account previously. In combination with the factorization formalism of non- 
relativistic QCD |l 8 ), this framework also leads to a successful description Q 
of inclusive J/if) and ?/>' production from S-meson decay JilJ. A rigorous pro¬ 
cedure to implement the finite-TOf, terms in an NLO framework where terms 
of the form a " +1 In" (/i 2 /m 2 ) are resummed by AP evolution and the univer¬ 
sality of the FFs is guaranteed by the factorization theorem [Q is to directly 
subtract the would-be collinear singularities of the form a s In (p|./m 2 ) in the 
massive-scheme result in a way that conforms with the massless scheme [^ 2 | 
(see also fll|). 
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